Connecting flow structures and heat flux in turbulent Rayleigh-Benard convection 
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The aspect ratio (T) dependence of the heat transfer (Nusselt number Nu in dimensionless form) 
in turbulent (two-dimensional) Rayleigh-Benard convection is numerically studied in the regime 
0.4 < T < 1.25 for Rayleigh numbers 10 7 < Ra < Ra 9 and Prandtl numbers Pr = 0.7 (gas) and 4.3 
(water). Nu(F) shows a very rich structure with sudden jumps and sharp transitions. We connect 
these structures to the way the flow organizes itself in the sample and explain why the aspect ratio 
dependence of Nu is more pronounced for small Pr. Even for fixed T different turbulent states 
(with different resulting Nu) can exist, between which the flow can or cannot switch. In the latter 
case the heat transfer thus depends on the initial conditions. 



PACS numbers: 47.27.-i, 47.27.te 

The physicists' view on fully developed turbulence has 
been dominated by Kolmogorov's seminal work [1 , pos- 
tulating universality of the small scales. Real- world flows 
however have walls and boundaries. How do these geo- 
metric properties affect the global transport characteris- 
tics of the flow such as heat or momentum transfer? The 
common view is that in the fully turbulent state, due 
to the large fluctuations of the system, the phase space 
is fully explored by the dynamics of the flow and that 
universality implies an at most weak dependence on the 
boundary conditions, clearly without any jumps in global 
transport properties. 

Only recently the community became aware of the pos- 
sibility of the coexistence of different turbulent states, 
with first or second order phase transitions in between 
them. Examples include different turbulent states in von 
Karman flow ("French washing machine" [2 ), in mag- 
netohydrodynamic turbulence in a closed system (von 
Karman turbulent liquid sodium experiment) [3], in tur- 
bulent rotating Rayleigh-Benard (RB) flow [4 , or in 
turbulent rotating spherical Couette flow [5], and pos- 
sibly even in high Rayleigh number turbulent RB flow 
in a cylindrical cell of aspect ratio T = 1/2 [6HT2] and 
T = 0.23 [13 . That transitions between different flow 
states occur in closed flows at low degree of turbulence 
has been well known and explored in the context of 
spatial-temporal chaos and pattern formation, see e.g. 
[14], but that they also occur between turbulent states 
clearly came as a surprise. 

In this paper we will numerically show that different 
turbulent states with first and second order phase tran- 
sitions in between them and even coexistence of different 
turbulent states also exist in an even simpler geometry, 
namely in a two-dimensional (2D) RB sample. The ad- 
vantages of the 2D geometry are (i) that a flow visualiza- 
tion is much easier, (ii) that the complicated axial and 
torroidal dynamics of 3D RB flow [16 does not compli- 
cate or even obscure the flow field analysis, and (iii) that 
it is numerically cheaper so that a good resolution in the 



aspect ratio T (i.e. many different runs with different T) 
becomes feasible. These features will enable us to reveal 
the connection between the flow organization and the 
heat-transfer properties. However, we stress that these 
connections also exist in 3D RB flow, and working them 
out from experimental measurements had been pioneered 
by Xia and coworkers [T5] . 

The main control parameter of this study will be the 
aspect ratio T, whose effect on the turbulent RB flow 
has not yet been extensively studied, with only a few 
theoretical and numerical studies [T7HT9] . The main rea- 
son for this short-coming is that it is experimentally very 
difficult to vary T in 3D cylindrical samples, as each T 
requires a new sample. That is why in experiment the T 
resolution has hitherto been insufficient to detect transi- 
tions between different turbulent states. In their review 
Ahlers et al. [16] conclude, based on a small number of 
experiments with few different Y [20-22], that the "weak 
T-dependence [of the heat transfer] suggests an insensi- 
tivity to the nature of the large-scale convection". We 
will show that in 2D RB this is definitely not the case. 

One may argue that the flow phenomena in 3D are 
different and richer, which is correct. However, various 
2D RB flow simulations [23 - 30 have revealed the value of 
this approach for a better understanding of turbulent RB 
convection and the analogies to 3D RB flow. In ref. [29] 
we could even reveal a one-to-one analogy between the 
flow organization in our 2D numerical simulations and 
in experiments in a quasi-2D sample of the same aspect 
ratio T = 1. 

The code on which the results in this paper are based 
is a fourth order finite difference discretization of the 
incompressible Oberbeck-Boussinesq equations. It has 
been described and tested in detail in ref. [27 . The ve- 
locity boundary conditions on the walls are no-slip. The 
(relative) temperature is fixed at A/2 at the bottom plate 
and —A/2 at the top plate, with adiabatic side- wall con- 
ditions. The grid resolution obeys the strict criteria for- 
mulated in ref. [31] . The imposed initial condition is a 
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FIG. 1. (color online) Nu(T) for Ra = 10 8 and (a) Pr = 4.3 
and (c) Pr = 0.7 (averaging time > 10 3 £#). The ac- 
companying figures (b) and (d) show the respective relative 
time the system is in the single roll state (SRS), double-roll 
state (DRS), triple-roll state (TRS), and quadruple-roll state 
(QRS). Additionally, the extend of time that no roll state 
could be identified is indicated. 

single roll state with a linear temperature profile. The 
sample aspect ratio is T — D/H, where D is the sample 
width and H the sample height. Times are given in mul- 
tiples of large eddy turnover times := 47t/(|cj c (£)|), 
where uj c is the center vorticity. Next to T, the other 
control parameters are the Rayleigh number Ra (the di- 
mensionless temperature difference between top and bot- 
tom) and the Prandtl number Pr = v j k, where v and 
k are kinematic viscosity and thermal diffusivity, respec- 
tively. The system responds with the heat transfer from 
bottom to top (in dimensionless form: Nusselt number 
Nu) and the turbulence intensity (in dimensionless form: 
Reynolds number Re). 

In fig. [l] we show Nu(T) for Ra = 10 8 and Pr = 0.7 
and 4.3, resulting from a very large averaging time of 
more than 10 3 t# , for which we have checked statistical 
convergence. For the lower Pr the curve shows varia- 
tions in Nu of up to 35%, with abrupt jumps and sharp 
transitions. From visual inspection of flow movies we de- 
duce that these jumps are associated with transitions in 
the overall flow organization from one to two vertically 
stacked rolls (at T w 0.9) and from two to three rolls (at 
T « 0.55). We automize this flow analysis by measuring 
the zeros in the vorticity along the center vertical axis, 
i.e., the zeros in uo(x = D/2,z). A state is not counted 
if its life time is shorter than The results for the 
relative time in a certain state is shown in fig. [I] For 
Pr = 0.7 this analysis confirms that the transitions are 
rather sharp, with the coexistence of different roll states 
only in a small T-range. Next to the jumps, relatively 
sharp transitions in Nu(T) can be observed. As an ex- 
ample we discuss the 35% increase in Nu from T = 0.45 
to T = 0.55, namely from Nu = 16 to Nu = 24. This 




FIG. 2. (color online) Snapshots of the temperature field for 
Ra = 10 8 and Pr = 0.7 for (a) V = 0.45 with the flow in 
a very stable TRS and (b,c) F = 0.55 with the flow in an 
unstable TRS. The top and bottom rolls in (b) grow until 
collision (c), which causes enhanced mixing and heat trans- 
port. Corresponding videos can be seen as supplementary 
material. The positions of the roll centers are automatically 
tracked and their PDFs as function of (x,z) are shown in (d) 
for T = 0.45 and (e) for T = 0.55. 

coincides with a developing instability of the triple roll 
state (TRS): The sample becomes too wide for three rolls 
so that the three-roll structure from time to time breaks 
down. This breakdown shuffles warm fluid upwards and 
cold fluid downwards, leading to a larger temperature 
gradient across the thermal boundary layers (BLs) and 
thus to increased heat flux. Typical snapshots of the 
temperature field in this regime are shown in fig. [2] the 
movies can be seen in the supplementary material. 

The enhanced vortex mobility at T = 0.55 as compared 
to T = 0.45 can also be deduced from the positions of 
the vortex centers, which we identify by a non-iterative 
ellipse-fit [32] coupled with a vortex criterium based on 
the non-zero imaginary part of the velocity gradient ten- 
sor's eigenvalues. The 2D probability functions (PDFs) 
of these positions are visualized in fig. |2]i (for T — 0.45) 
and fig. [2^ (for T = 0.55). While in the former case we 
see three pronounced sharp peaks, reflecting the stability 
of that flow configuration, in the latter case the peaks 
are smeared out, showing the vortex mobility which con- 
tributes to the enhanced heat flux. 

We now come to the case of large Pr = 4.3. Here the 
transitions in Nu(T) are much smoother as compared to 
the Pr = 0.7 case (compare figs. [TJd and d) with wide T 



FIG. 3. (color online) Two temperature snapshots for Ra = 
10 9 , Pr = 4.3, and T = 0.93. The system is either in the SRS 
(a) or DRS (b). The corresponding movie is available online 
as supplementary material. Heat transport efficiency differs 
between these states, see fig. [4] 

regimes of coexistence between different states, between 
which the dynamical evolution of the flow meanders, as 
e.g. seen from temperature field snapshots of fig. [3j the 
supplementary movie, and the time series of the Nus- 
selt number. For these larger Pr the dynamics is much 
richer as the flow keeps on switching between the different 
states. The reason for this lies in the plume dynamics: 
Whereas for Pr = 0.7 the larger thermal diffusivity tends 
to smoothen out temperature differences between differ- 
ent flow parcels, for larger Pr = 4.3 the plumes keep 
their thermal energy for a longer time. Thus the plumes, 
which feed the competing rolls as described in ref. [29] , 
have a longer life time, which leads to more lively roll 
dynamics. Therefore a 100% dominance of a certain roll 
state, common for small Pr, does not occur. 

Based on the state analysis, the heat flux can be condi- 
tionally averaged on the SRS, DRS, or TRS. The results 
for these conditional heat fluxes are shown in fig. [4] For 
the vertically arranging vortices it again holds that the 
stable states with n rolls enable larger heat transfer than 
the stable states with n+1 rolls. For n=l and n=2 the 
difference is about 5% for Ra = 10 9 , Pr = 4.3. Re- 
markably, that value is similar to what one would expect 
for the 3D situation, based on an extrapolation towards 
Ra = 10 9 of the Nu measurements for the SRS and DRS 
in ref. [33] (for 4 • 10 9 < Ra < 10 11 ), which roughly gives 
3%. 

We also found a case for which the final state and thus 
Nu depend on the initial conditions of the flow field, see 
fig. [H] for Ra = 10 7 and Pr = 0.7, which according to 
ref. [25] is already turbulent. The unstable TRS initial 
condition for T = 0.63 falls back to the DRS, whereas 
for initial DRS condition at T = 0.64 the flow remains 
in this DRS (with lower Nu) for longer than 3000 tg?. It 
remains to be seen whether this explicit dependence on 
initial conditions can also be found at larger Ra. 

We finally address the question why Nu is much more 
sensitive to the flow organization at small Pr as com- 



FIG. 4. (color online) Conditionally averaged Nu(T) for the 
various states SRS, DRS, and TRS, which can be identified 
at Ra = 10 9 and Pr = 4.3. For fixed T the flow jumps in 
between these different states. 
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FIG. 5. (color online) Nu(T) for Ra = 10 7 and Pr = 0.7. 
For r « 0.64 the Nusselt number depends on the initial con- 
ditions, see text for details. 

pared to large Pr, see again fig. [I] The first reason for 
this finding is the organization of the thermal BL and 
kinetic BL: For large Pr the thermal BL is nested in 
the kinetic one. Thus modifications of the bulk flow are 
buffered by the kinetic BL and hardly lead to a differ- 
ent thermal BL thickness \q and therefore to hardly any 
change in Nu « L/(2\q). In contrast, for low Pr the 
thermal BL is thicker than the kinetic one and is there- 
fore fully exposed to bulk flow modifications, leading to 
a much stronger dependence of Nu ~ L/(2\q) on the 
flow organization. The second reason for the larger sen- 
sitivity of Nu at smaller Pr is the larger thermal diffu- 
sivity, leading to a thermal smoothening of the plumes. 
As discussed above this results in sharper transitions be- 
tween different states. Our interpretation is supported by 
fig. |6j in which we plot Nu(T) vs a normalized Re z (T), 
where Re z = u z ^ rrns L/v ', which represents the bulk flow. 
For Pr = 4.3 there is hardly any dependence of Nu on 
Re z j (Re z ) , whereas for Pr = 0.7 this dependence is ma- 
jor. Here (Re z ) is the average Reynolds number of the 
shown T range 0.4 < T < 1.25. 
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FIG. 6. (color online) Nu(T) vs. Re z (Y)/(Re z ) for Pr 
(red squares) and Pr = 4.3 (blue bullets) at Ra — 10 8 . 
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In conclusion, we could clearly identify different turbu- 
lent states, corresponding to different roll structures. It 
is remarkable that these features, which have been associ- 
ated to lower Reynolds and Rayleigh number flow in the 
regime of spatial-temporal chaos with much less dynami- 
cal degrees of freedom, survive for such high Ra. It seems 
that the turbulent states 'live' on low-dimensional struc- 
tures. The different turbulent states are associated with 
different overall heat transfer. Transitions between the 
states - either dynamically or as function of the control 
parameter V - thus imply jumps and sharp transitions in 
Nu(T). Based on the flow organization, we understand 
why these variations are larger for lower Pr and why the 
overall flow dynamics is richer at larger Pr, for which 
the plumes can keep their thermal identity for a longer 
time. The next step will be to push the present results 
to much higher Ra in order to see whether these features 
survive. Based on our analysis for Ra = 10 8 and 10 9 we 
presume that this could be the case as the stability of the 
states and the variation in Nu(T) is higher for Ra = 10 9 
compared to Ra = 10 8 at Pr = 4.3. Additionally, large 
scale roll states have been distinguished up to Ra = 10 10 
[29] and can probably be found for even higher Ra. 
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